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We study the electronic properties of a two-dimensional quasiperiodic tiling, the isometric general- 
ized Rauzy tiling, embedded in a magnetic field. Its energy spectrum is computed in a tight-binding 
approach by means of the recursion method. Then, we study the quantum dynamics of wave packets 
and discuss the influence of the magnetic field on the diffusion and spectral exponents. Finally, we 
consider a quasiperiodic superconducting wire network with the same geometry and we determine 
the critical temperature as a function of the magnetic field. 
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I. INTRODUCTION 



For periodic systems such as the square, the triangular 
or the hexagonal lattice embedded in a magnetic field, 
it is well-known that commensurability effects between 
the characteristic magnetic length and the lattice period 
may lead to very interesting features. For example, if 
we consider a standard tight-binding model, the energy 
spectrum displays a complex multifractal structure when 
varying the magnetic flux per unit cell as illustrated by 
the famous Hofstadter butterfly ||, ||]. However, this 
spectrum is strongly geometry-dependent since even for 
periodic tilings, some anomalies can appear leading to 
interesting quantum interference phenomena |Q, ||. This 
dependence has led several groups to analyze the inter- 
play between the quasiperiodicity and the magnetic frus- 
tration 1^, 0, ^. Nevertheless, in these studies, it is 
difficult to know the real contribution of the quasiperi- 
odic order since the incommensurability between the tiles 
area is another source of aperiodicity in the spectrum Q . 
We must indeed distinguish several classes of tilings de- 
pending on their structural order (periodic, quasiperiodic 
or disordered), their topology (trivial or nontrivial) and 
the ratio of the tiles area (rational or irrational). If we 
focus on quasiperiodically ordered structures, there is ac- 
tually one type of tiling (those of codimension one) that 
has never been studied so far which consists in a non- 
trivial topology (several types of coordination numbers) 
with commensurate tile areas. 

The aim of this work is to study such tilings by 
considering a modified version of the two-dimensional 
generalized Rauzy tilings (GRT) presented in Ref. 
that are built with tiles having the same areas. Note 
that a great advantage of such systems is that all the 
physical quantities are periodic functions of the magnetic 
field so that we can restrict our analysis to a finite range 
of flux values. 



This paper is organized as follows: in the next sec- 
tion, we briefly present the two-dimensional isometric 
generalized Rauzy tilings (iGRT) whose properties will 
always be compared to the square lattice ones. In Sect. 
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III, we introduce the tight-binding Hamiltonian and we 
compute the energy spectrum of both tilings. Wc show 
that despite the quasiperiodic order, the spectrum of the 
iGRT has a fine nontrivial structure. To get more pre- 
cise informations on the spectrum and the eigenstates, 
we compute in Sect. EV the time evolution of wave pack- 
ets for different fluxes per plaquette 4> and we focus on 
the average autocorrelation function and the mean square 
spreading. For the iGRT a singular continuous spectrum 
and a sub-ballistic propagation is found for rational and 
irrational / = 4>/4>{) where 4)q is the flux quantum. By 
contrast, for the square lattice, the propagation is bal- 
listic and the spectrum absolutely continuous for ratio- 
nal /, whereas the spectrum is singular continuous and 
the spreading sub-ballistic for irrational /. Sect. ^ is 
devoted to the superconducting wire network with the 
iGRT geometry. After mapping the linearized Ginzburg- 
Landau equations onto a tight-binding model with non- 
constant hopping terms, we compute the transition tem- 
perature as a function of the magnetic field and compare 
it with those of the square lattice. A cusp-like struc- 
ture is clearly observed and shows the importance of the 
quasiperiodic order in this context. 



II. THE ISOMETRIC GENERALIZED RAUZY 
TILINGS 

The GRTs are codimension one quasiperiodic struc- 
tures that can easily be built in any dimension D using 
the standard cut and project method These tilings 
have a complex topological structure with sites of coor- 
dination number ranging from Z) -I- 1 to 2D -f 1 . In this 
study, we focus on the two-dimensional GRT whose con- 
truction is based on the irrational solution 9 ~ 1.839 of 
the equation: 

x^^^x^ + x + l. (1) 

Contrary to previous studies |l^, ^ concerning these 
tilings, we consider here an isometric version of the 2D 
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GRT where all the edge lengths are equal. This can be 
very easily done by projecting the selected sites of the 
3-D cubic lattice perpendicularly to the direction (1,1,1) 
instead of the Rauzy direction {1,9,6^). Such a choice 
is motivated by the fact that since these tilings are codi- 
mension one tilings, equal edge lengths imply equal tile 
areas. This property which is only true for codimension 
one tilings allows us to concentrate on the role played 
by the topological quasiperiodicity related to the connec- 
tivity matrix. 




FIG. 1: A piece of the 2D isometric generalized Rauzy tiling. 
All the tiles have the same area. 



III. HAMILTONIAN AND BUTTERFLIES 

In this study, we consider a standard tight-binding 
Hamiltonian given by: 

H = -J2u,\i){jl (2) 

where |i) is a localized orbital on site i. The hopping term 
tij = 1 if i and j are nearest neighbors and otherwise. 
In the presence of a magnetic field [Q , tij is multiplied 
by a phase factor e*"*"'^ involving the vector potential A: 

^ = / A.dl, (3) 

where (j)Q — hc/e is the flux quantum. In the follow- 
ing, we only consider the case of a uniform magnetic 
field B = Bz which can be obtained, for example, with 
the Landau gauge A = i3(0,a:,0), and we denote by 
(j) ~ B(j?\f%j2 the magnetic flux through an elementary 
tile. It is also useful to introduce the dimensionless pa- 
rameter / = ^/(/>o known as the reduced flux. 

The study of the zero-field spectrum of H has been 
presented in Ref. |ll|, |l^, ^ for the nonisometric GRTs, 
but for _B = the tile areas is an irrelevant parameter 



since H is, in this case, simply proportionnal to the con- 
nectivity matrix. The spectral measure is absolutely con- 
tinuous, the spectrum is gapless, and the weakly critical 
eigenstates, are responsible for an anomalous diffusion. 
To obtain the spectrum of R for i? 7^ 0, wc have used 
the recursion method also known as Lanczos algo- 
rithm ||l^ . This powerful tool allows one to compute the 
local density of states (LDOS) of a given initial state. 
The main advantage of this method is that we can eas- 
ily compute this LDOS for very large system size (about 
10^ sites), by diagonalizing only a small half-chain of a 
few hundred sites. However, it does just give a local in- 
formation but, if one considers an initial random phase 
state, its LDOS gives a fairly good approximation of the 
total density of states. To compute these spectra we have 
computed 500 recursion coefficients in a circular cluster 
of 785757 sites for the square lattice, and 460 coefficients 
in a circular cluster of 767604 sites for the iGRT. We 
must also mention an important point concerning the 
edge states. Indeed, since we compute the recursion co- 
efficients up to a given order, there are obviously some 
spurious states that arise due to the fact that the recur- 
sion cluster is finite with open boundary conditions. To 
get rid of them, we have selected eigenvectors for which 
the amplitude on the initial state of the recursion was 
greater than a given threshold (here 10~^). 



As it can be seen in Fig. ||, the spectrum of the square 
lattice obtained here with the recursion scheme is in very 
good agreement with those computed by Hofstadter in 
1976 using the mappping onto the Harper equation [0. 
As it can be readily seen, both spectra display the sym- 
metry E ^ —E related to the fact that both tilings are 
bipartite. In addition, since there is only one characteris- 
tic area, the Hamiltonian is a periodic function of / with 
period 1. Finally, the symmetry with respect to / = 1/2 
just reflects the symmetry tij ~tij- At small fields, 
the emergence of Landau levels, expected in the contin- 
uum limit, is clearly observed. In the square lattice, it 
is well known that rational values of / = p/q lead to an 
absolutely continuous spectrum made up of q subbands 
that can touch each other as for /= 1/2. Of course, 
nothing similar is expected for quasiperiodic tilings since 
no periodicity exists. Nevertheless, it is interesting to note 
that for the iGRT, there is still a complex gap structure 
with notably a large gap that can be followed continu- 
ously from / = (where it vanishes) till / = 1/2. We 
can also see that the iGRT spectrum display a fine sub- 
band structure that would require a detailed but difficult 
analyzis since even in the zero field case, it is nontrivial. 



To go beyond this basic description, we have analyzed 
the nature of the spectrum for particular values of the 
reduced fiux by studying the quantum dynamics of wave 
packets in both structures for remarkable values of the 
reduced fiux. 
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FIG. 2: Energy spectrum as function of the reduced flux / = 
(j)/4>o for the square lattice (top) and for the iGRT (bottom). 
For convenience, we have only considered rational values of 
/ = p/q for q < 37. 



IV. WAVE PACKETS DYNAMICS 

The quantum dynamics in two-dimensional quasiperi- 
odic tilings has been widely investigated within the last 
ten years either by exact diagonalizations p^, [l9| o r 
by approximating the evolution operator e^^* ||ll[7l2[. 



Here we have used the Second Order Differencing Scheme 
pO| , that consists, for a given initial state \^{0)), to make 
the following approximation: 



\il;{t + At)) = \i;{t ~ At)) - 2iAtH\iP{t)). 



(4) 



Contrary to exact diagonalizations, this method allows 
one to investigate the dynamics on very large systems 
(about one million sites). Practically, choosing a time 
step At — 0.05 is sulBcient to get a good accuracy and 
all our results have been obtained with this value. We 
have considered two quantities of interest: the averaged 
autocorrelation function 



C{t) ^-j dt'P{t') - 1- 
t Jo 



(5) 



where P{t) = \{tlj{0)\ip{t))\'^ and the mean square spread- 
ing 



where R is the position operator, the origin being taken 
at the center of the cluster. The spectral and diffusion ex- 
ponents a and (3 are defined from the long time behavior 
of C and AR^ respectively. If a = I, the spectrum is ab- 
solutely continuous, if a = 0, the spectrum is pure point, 
and if < a < 1 the spectrum is singular continuous. 
Concerning the diffusion exponent, [3—1 corresponds 
to a ballistic propagation whereas < /5 < 1 defines a 
sub-ballistic regime {(3 = 1/2 is the diffusive case). 

For the square lattice, since all sites are equivalent, 
we have considered an initial state localized on a sin- 
gle (central) site. For the iGRT, the situation is more 
complicated since, as for the Octagonal tiling |l^, the 
exponents depends, a priori on the initial conditions as 
shown in Ref. |ll|, Of course, this is only true at 

short times, since in the asymptotic regime, we expect 
a = inf£;{a(i?)}, /3 = sup^{/3(£')} for E in the spec- 
trum. Here, we have built an initial random phase state 
localized on a circular cluster 10^ sites) whose radius 
is about 1/100 of the cluster considered for the compu- 
tation, and we have checked that our results were weakly 
sensitive to the precise random phase configuration cho- 
sen. However, we can never ensure that we have reached 
the asymptotic regime so that the exponents given here 
just give a qualitative idea of the wave packet spreadings. 

We have considered several special values of the re- 
duced flux: (i) / = for which the band width are max- 
imum for both lattices ; (ii) J — 1/2 for which the band 
width of the iGRT is minimum ; (iii) / = for which 
the band width of the square lattice is minimum ; (iv) 
f — which is a natural irrational number in this con- 
text. Our results are plotted in FigJ| and Fig.^. 

For the square lattice, everything is qualitatively 
known for these values. Indeed, when / is a rational 
number, the eigenstates are Bloch waves, the spectrum 
is absolutely continuous (a = 1) and the spreading is 
ballistic {(3 = 1). For f — 0, when starting from a single 
site, one can even easily show that: P{t) = Jq (2t) which 
behaves as t~^ at large time so that C{t) '--^ t^^ and 
AR^{t) = At^. For irrational /, the spectrum is known 
to be singular continuous |l8, 21 , and the propagation 
sub-ballistic Nevertherless, we emphasize that 

existing results about the dynamics have been obtained 
by analyzing the Harper chain and do not reflect the dy- 
namics on the square lattice. Indeed, the mapping of the 
square lattice onto the Harper chain results from a de- 

■4;{x) 
Thus, 
I only 



AR'it) = {m\R'\m)~mt)\R\m)? 



2/3 



(6) 



coupling of the eigenstates given by ^{x, y) 
which is due the special Landau gauge choice, 
the dynamics discussed in Refs. ||l^, H^, |_ 
concerns the x direction (at fixed ky) and thus do not 
give the exponent of the two-dimensional square lat- 
tice. For example, we have computed the dynamics for 
/ = (1 + VZ)/2 = T and we have found a — 0.50(1) in 
the square lattice instead of a 0.14(1) in the Harper 
chain ||, m. 

Interestingly, for / = \/2, we have found a diffu- 
sive propagation with (3 = 0.50(1) whereas for f = 6, 
the problem is more complicated since there are clearly 
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FIG. 3: Averaged autocorrelation function. Top: square 
lattice with / = 0, 1/2, ^p2,Q (from bottom to top). Bottom: 
iGRT for / = 0, 1/2, e, \/2 (from bottom to top). 



FIG. 4: Mean square spreading. Top: square lattice for 

/ = 0, 1/2,6*, (from top to bottom). Bottom: iGRT for 

/ = 0, 1/2, 6*, \/2 (from top to bottom). 



two regimes. If we consider only t > 250 we obtain, 
/3 — 0.94(2) (weakly sub-ballistic motion) but we can- 
not decide whether it is an intermediate regime or not. 
The exponents obtained by fitting our data are given in 
Table I. 

For the iGRT, the situation is different since even for 
f — 0, the spectrum has some critical states responsi- 
ble for a sub-ballistic spreading jl^. However, as shown 
in Ref. the diffusion exponent is strongly energy- 
dependent and might even be very close to 1 at the 
band edges. In the asymtotic regime, a quasi-ballistic 
propagation should therefore be recovered. Contrary to 
the square lattice, a sub-ballistic spreading is found for 
/ = 1/2 with f3 = 0.61(1) which is not really surprising 
since even though 1/2 is a rational number, the underly- 
ing quasiperiodicity of the iGRT does not allow for the 
contruction extended eigenstates. In a sense, there are 
no commensurability effects for rational /. For / = a/2, 
it is interesting to note that the propagation is faster 
in the iGRT {j3 = 0.57(2)) than in the square lattice 
f3 = 0.50(1). This means that the quasiperiodic order 
can break the destructive quantum interference effects 
that slow the spreading down for incommensurate / in 
periodic structures. Nevertheless, this is not always the 
case, since for f = 9, the propagation is faster in the 
square lattice than in the iGRT. Concerning the nature 
of the spectrum, it is always found to be clearly singular 
continuous but for / = 0, we cannot be sure that a ^ 1 



in the asymptotic regime. 

V. SUPERCONDUCTING WIRE NETWORKS 

In this section, we analyze the role of the quasiperiodic 
order on the superconducting-normal transition temper- 
ature as a function of the magnetic field Tc{f). As ex- 
plained above, the main advantage of the iGRT is that 
all its elementary tiles have the same area. Thus, we can 
really analyze the role of the topological quasiperiodic 
order which has, in fact, never been investigated. 

The study of quasiperiodic wire networks has been ini- 
tially motivated by several experiments |2j, |2^, |2^, |2^ 
in which the superconducting transition temperature as 
a function of the magnetic field was determined. Using a 
mapping between the linearized Ginzburg-Landau equa- 
tions and the tight-binding model, the rich structure of 
this phase boundaries has been then computed by Nori 
and coworkers j2^, |2^ . The main idea of this correspon- 
dence developped by Alexander and De Gennes |3|] 
is to integrate the LGL equations on each wire and to 
obtain a set of equations that determines the order pa- 
rameter ij^i at each nodes i. This can be simply achieved 
using the continuity of the order parameter and the cur- 
rent conservation at each node. Denoting by Uj the dis- 
tance between the (connected) nodes i and j and by ^ the 
coherence length, the current conservation at each node 
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0.57(2) 
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TABLE I: Spectral and diffusion exponents for the square 
lattice and for the iGRT. 



I expresses as: 



Tpl ^ COt(/j 



'/sin(;,,/O=0, (7) 



where the elementary flux quantum involved on the defi- 
nition (H) of 7y is here 4>'q = hc/2e. The coherence length 
depends on the temperature: 



TcjO) 



1/2 



(8) 



where Tc{0) is the transition temperature at zero field. 
For a given value of the magnetic field, the transition 
temperature is thus given by the highest temperature 
for which a nontrivial solution of exists. If we now 
consider a system where all the lengths are equal kj = I, 
and if we set (j^i = y/Zitpi where Zi is the coordination of 
the node i, Eq. (0) can be rewritten as ||29|| : 



cos{l/^)(j)i 



(9) 



Setting with tij — —{ziZj)^^/^ and E = cos(//^), we ob- 
tain the secular system given by the Hamiltonian (H) and 
thus the field dependence of the transition temperature: 



1 



Tcif) 

TciO) 



(10) 



where Eo{f) is the highest eigenvalue of this modified 
Hamiltonian for a reduced flux / = 27r(/()/0Q. 

For the square lattice, the first experiments have been 
performed by Pannetier et al. and have confirmed 
the validity of these predictions. Indeed, a self-similar 
cusp-like structure for rational / has been clearly iden- 
tified in agreement with the fractal pattern of Eo{f) al- 
ready pointed by Hofstadter For quasiperiodic sys- 
tems, experiments by Behrooz and coworkers |^, |25[] 
have also revealed an cusp-like structure that can be 
understood in terms of incommensurability between the 
tiles area 
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We have computed Tc(/) for the square lattice and for 
the iGRT using, as previously, the recursion method. As 
explained in Ref . , this is a particularly well-adapted 
algorithm to obtain the ground state of the modified 
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FIG. 5: Reduced critical temperature as a function of the 
reduced flux / for the square lattice (grey curve) and for the 
iGRT (black curve). Inset: zoom around f — 1/2 for the 
iGRT. 



Hamiltonian which is in this case equal to —Eo{f). The 
results displayed in Fig. ||, clearly show the strong cusp- 
like structure for the square lattice but nothing similar 
seems to occur for iGRT. However, when zooming, we 
can distinguish a fine structure in Td f) but with oscilla- 
tions much smaller than those of the square lattice. We 
also note that, for the iGRT, the minimum of Tc{f) is ob- 
tained for f — 1/2 where an inverse-cusp is observed as 
for the Dice lattice 1 33 or for the Kagome lattice ^ . 
This rich cusp structure is rather strange for this tiling 
since there is no incommensurability effect due to either 
the ratio of the tile areas (which are all equal for the 
iGRT) or to the lattice periodicity. 



VI. CONCLUSION 

We have studied the spectrum and the quantum dy- 
namics of an electron in a quasiperiodic tiling embed- 
ded in a magnetic field. As in periodic structures, its 
spectrum displays a complex energy pattern when the 
magnetic flux is varied. We have shown that for several 
values of the reduced flux /, the wave packet spreading 
is sub-ballistic with a field-dependent diffusion exponent. 
For / = we have even found a propagation which is 
faster in the quasiperiodic tiling than in the square lat- 
tice for which the motion is, in this case, diffusive. It 
would be very interesting to study the influence of the 
disorder (topological, geometrical, or Anderson-like) on 
these exponents since, even in the zero-field case, non- 
trivial behaviours have already been observed p^ . 

We have also determined the superconducting-normal 
transition line Tc{f) that displays a cusp- like structure 
revealing the importance of the quasiperiodic order. 
However, the analysis of the cusp positions clearly 
requires further investigations. It would be also im- 
portant to study the transition line for codimension 
one random tilings that also have one type of tile but 
no quasiperiodic order. In this case, the absence of 
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cusp would confirm the particular role played by the 
quasiperiodicity. 
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